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Abstract. We have demonstrated that the isospin of nuclei influences the fragment production during the 
nuclear liquid-gas phase transition. Calculations for Au 197 , Sn 124 , La 124 and Kr 7S at various excitation 
energies were carried out on the basis of the statistical multifragmentation model (SMM). We analyzed 
the behavior of the critical exponent r with the excitation energy and its dependence on the critical tem- 
perature. Relative yields of fragments were classified with respect to the mass number of the fragments in 
the transition region. In this way, we have demonstrated that nuclear multifragmentation exhibits a 'bi- 
modality' behavior. We have also shown that the symmetry energy has a small influence on fragment mass 
distribution, however, its effect is more pronounced in the isotope distributions of produced fragments. 

PACS. 25.70.Pq Multifragment emission and correlations - 21.65.+f Nuclear matter - 24.60.-k Statistical 
theory and fluctuations 



1 Introduction 

One of the most interesting phenomena in heavy ion re- 
actions is multifragmentation of nuclei, which is a very 
promising process to determine the properties of nuclear 
matter at subnuclear densities and high excitation ener- 
gies. Experimental and theoretical analysis of this phe- 
nomenon is also very important for our understanding of 
astrophysical events, such as supernova explosions and for- 
mation of neutron stars [T|P]. The ground state properties 
of nuclear matter can usually be determined as a theo- 
retical extrapolation, based on nuclear models designed 
to describe the structure of real nuclei. When nuclei are 
excited, these ground state properties show small changes 
without structural effects. At very high excitation energies 
(E* > 2 — 3 MeV/nucleon) the nuclei can slowly expand 
by remaining in a state of the thermodynamic equilibrium. 
During this expansion, they enter the region of subsatu- 
ration densities, where they become unstable to density 
fluctuations, and break up into the fragments (droplet for- 
mation). In this case, it is believed that a liquid-gas type 
phase transition is manifested. Even though there exist 
some theoretical and experimental evidences of this phase 
transition, so far some uncertainties remain in its nature 
PHH|3| |5]17ll51l^lllUllllj . Recently, the problem of isospin ef- 
fects in multifragmentation gains a lot of interest in con- 
nection with the astrophysical applications. The main goal 
of this paper is to analyze these effects in the framework 
of a realistic model. 



2 Statistical multifragmentation model 

A model used for extracting the properties of nuclear mat- 
ter, should be capable of reproducing the observed frag- 
mentation properties. Statistical multifragmentation model 
(SMM) |12U13| has provided a good reproduction of exper- 
imental data, as shown by many experimental analyses 
9 13 14. 15. l 6Tll7| . Within the SMM we consider the whole 
ensemble of breakup channels consisting of hot fragments 
and nucleons in the freeze-out volume. In finite nuclear 
systems the SMM takes into account the conservation laws 
of energy E* , momentum, angular momentum, mass num- 
ber Aq and charge number Z$. An advantage of SMM is 
that it considers all break-up channels including the com- 
pound nucleus and one can study the competition among 
them. In this respect, there is a natural connection be- 
tween multifragmentation and traditional decay channels 
of nuclei as the evaporation and fission at low excitation 
energies. On the other hand the SMM can also address the 
thermodynamical limit, i.e. very large systems, that makes 
it suitable for astrophysical applications ■ The probabil- 
ity of any breakup channel is assumed to be proportional 
to its statistical weight, and the channels are generated 
by the Monte Carlo method. The physical quantities such 
as average mass and charge fragment yields, temperature 
and its variance are calculated by running the summations 
over all members of the ensemble. 

In this work we use the SMM version described fully 
in |13| . However, for this paper, it is important to ex- 
plain how the hot primary fragments are treated. Light 
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fragments with mass number A < 4 and charge num- 
ber Z < 2 are considered as stable particles ("nuclear 
gas") with masses and spins taken from the nuclear ta- 
bles. Only translational degrees of freedom of these par- 
ticles contribute to the entropy of the system. Fragments 
with A > 4 are treated as heated nuclear liquid drops, and 
their individual free energies F A z are parameterized as a 
sum of the bulk, surface, Coulomb and symmetry energy 
contributions 

F AZ = F B AZ + F S AZ + E C AZ + E%£. (1) 

The bulk contribution is given by Ff z = (-W -T 2 /e )A, 
where T is the temperature, the parameter eo is related to 
the level density, and Wq = 16 MeV is the binding energy 
of infinite nuclear matter. Contribution of the surface en- 
ergy is F S AZ = B A 2 / 3 (g^) 5 / 4 , where B = 18 MeV 
is the surface coefficient, and T c = 18 MeV is the criti- 
cal temperature of the infinite nuclear matter. Coulomb 
energy contribution is E AZ = cZ 2 /A 1 / 3 , where c is the 
Coulomb parameter obtained in the Wigner-Seitz approx- 
imation, c = (3/5)(e 2 /r )(l — (p/po) 1 ^ 3 ), with the charge 
unit e, ro =1.1 7 fm, and po is the normal nuclear matter 
density ~ 0.15/to~ 3 . And finally, the symmetry term is 
E* A V £ = 7 (A - 2Z) 2 /A, where 7 = 25 MeV is the symme- 
try energy parameter. All the parameters given above are 
taken from the Bethe-Weizsacker formula and correspond 
to the assumption of isolated fragments with normal den- 
sity in the freeze-out configuration. This assumption has 
been found to be quite successful in many applications. 
However, for a more realistic treatment one should con- 
sider the expansion of the primary fragments in addition 
to their excitations. The residual interactions among them 
should be considered as well. These effects can be taken 
into account in the fragment free energies by changing the 
corresponding liquid-drop parameters. 

Here, we have applied the SMM model to nuclear sources 
of different mass and isospin which can easily be used in 
modern experiments. The excitation energy range is of 
E* = 2 — 20 MeV/nucleon, and the freeze-out volume, 
where the intermediate mass fragments are located after 
expansion of the system, is V = 3Vq (Vb = 47rylor 3 l /3 
is the volume of a nucleus in its ground state). In other 
words, the fragment formation is described at a low den- 
sity freeze-out (p w po/3), where the nuclear liquid and 
gas phases coexist. The SMM phase diagram has already 
been under intensive investigations (see e.g. |18j ) with the 
same parametrization of the surface tension of fragments 
versus the critical temperature T c for the nuclear liquid- 
gas phase transition in infinite matter. Therefore, it is pos- 
sible to extract T c from the fragmentation data ^l- The 
symmetry energy contributes to the masses of fragments 
as well. Taking this contribution into account we shall 
investigate the influence of the symmetry energy coeffi- 
cient 7 on fragment production in the multifragmentation 
of finite nuclei. In this paper we concentrate mainly on 
properties of hot fragments. However, their secondary de- 
excitation is important for final description of the data, 
which was included in SMM long ago, and its effects were 
already discussed somewhere else f [T^lll9ll2()] V Here we 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
A/Ao A/A 



Fig. 1. Relative yield of hot primary fragments versus A/Ao 
for Au 197 , Sn 124 , La 124 and Kr 78, at different excitation ener- 
gies 3, 4, 5 and 8 MeV/nucleon. 

demonstrate how the secondary deexcitation code can be 
modernized so that the change in the symmetry energy of 
nuclei during the evaporation is taken into account. 

3 Isospin dependence of fragment 
distributions 

For our calculations, we consider Au 197 , Sn 12i , La 124 and 
Kr 78 nuclei to see how isospin affects the multifragmenta- 
tion phenomena with their neutron-to- proton ratios (N/Z), 
which read 1.49, 1.48, 1.18 and 1.17 respectively. The two 
sources Au 197 and Sn 124 have nearly the same N/Z ratios 
1.49 and 1.48, respectively, but different mass numbers, 
whereas Sn 124 and La 124 have the same mass numbers, 
but very different N/Z ratios 1.48 and 1.18, respectively. 
To see the effect of the size of the nucleus in our study, we 
have also included Kr 78 , which has the lowest mass num- 
ber among the considered nuclei here, but its N/Z ratio 
1.17 is very close to that of La 124 . Relative yield of hot 
fragments produced after the break-up of Au 197 , Sn 124 , 
La 124 and Kr 78 nuclei as a function of A/Aq is given in 
Fig. 1 for E* = 3, 4, 5 and 8 MeV/nucleon. As was shown 
in previous analysis of experimental data (see e.g. [711^1 
H3U15p these results are fully consistent with experimen- 
tal observations. The effect of isospin of different sources 
in liquid-gas phase transition region can be clearly seen 
for different excitation energies. The mass distribution of 
fragments produced in the disintegration of various nu- 
clei evolves with the excitation energy. As can be seen in 
Fig. 1, for Au 197 and Sn 124 multifragmentation onset is 
about 5 MeV/nucleon and for La 124 and Kr 78 it is about 
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Fig. 2. N/Z ratio of hot primary fragments versus A for Au 197 , 
Sn 124 , La 124 and _ft"r 78 at different excitation energies 3, 4, 5 
and 8 MeV/nucleon. 
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Fig. 3. The values of critical exponents r (top panel) and 
tz (bottom panel), as a function of excitation energy E* for 
Au 197 , Sn 124 , La 124 and Kr 78 . 



4 MeV/nucleon at standard SMM parameters so that U- 
shape mass distributions disappear at these energies. We 
can also describe this evolution with the temperature (see 
the caloric curves below). At low temperatures (T < 5 
MeV), there is a U-shape distribution corresponding to 
partitions with few small fragments and one big residual 
fragment. At high temperatures (T > 6 MeV), the big 
fragments disappear and an exponential-like fall-off is ob- 
served. In the transition region (T ~ 5 — 6 MeV), how- 
ever, one observes a transition between these two regimes, 
which is rather smooth because of finiteness of the sys- 
tems. 

In Fig. 2 we demonstrate N/Z ratios of the hot frag- 
ments produced in the freeze-out volume of the same sys- 
tems. One can see that neutron content of hot fragments 
is only slightly smaller than that of the whole system since 
the most of the neutrons are still contained in fragments. 
It is also seen that the neutron richness of fragments is 
increasing with their mass number |21j . This is a general 
behavior for nuclear systems in equilibrium, where binding 
energy is influenced by interplay of the Coulomb and the 
symmetry energy. For some systems one can observe spe- 
cific isotopic effects such as the increasing of the neutron 
richness of intermediate mass fragments (with Z=3-20) 
with the excitation energy (compare the results for 3 and 
8 MeV/nucleon). This is because after removing heavy 
fragments by increasing excitation energy, their neutrons 
are accumulated into the smaller fragments |2*T] . 

Usually a power-law fitting is performed with, Y(A) oc 
A- T and Y(Z) cx Z~ Tz , where Y(A) and Y{Z) denote 
the multifragmentation mass and charge yield as a func- 
tion of fragment mass number A and charge number Z, 



respectively. The parameters r (for mass distribution) and 
tz (for charge distribution) can be considered in thermo- 
dynamical limit as critical exponents 16,22 . In calcula- 
tions, we consider the fragments in the range 6 < A < 40 
for mass and 3 < Z < 18 for charge number. The lighter 
fragments are considered as a nuclear gas. The obtained 
values of r and Tz at T c = 18 MeV as a function of E* are 
given in Fig. 3 for cold fragmentation (with secondary de- 
excitation). There is a minimum at about E* = 5.3, 5.2, 4 
and 4.2 MeV/nucleon for Au 197 , Sn 124 , La 124 and Kr 7s , 
respectively. The results for Au 197 and Sn 12A are very sim- 
ilar to each other and significantly different from those 
obtained for La 124 and Kr 7s . This means that interme- 
diate mass fragment (IMF) distributions approximately 
scale with the size of the sources, and they depend on the 
neutron-to-proton ratios of the sources. This is because 
the symmetry energy still dominates over the Coulomb in- 
teraction energy for these intermediate-size sources. One 
may also see from these results that the lower N/Z ra- 
tio leads to smaller r and Tz parameters with the flatter 
fragment distribution. A large N/Z ratio of the source fa- 
vors the production of big clusters since they have a large 
isospin. This is the reason for domination of partitions 
consisting of small IMFs with a big cluster in the transi- 
tion region (U-shape distribution) . It is also seen from the 
bottom of Fig. 3 that the values of tz are very close to the 
values of r since the neutron-to-proton ratio of produced 
IMFs exhibits a small change within their narrow charge 
range (see also [T5ll2*T] ). 

For completeness, we have calculated the average IMF 
multiplicity and its variance for the excitation energy range 
of 2 — 20 MeV/nucleon for all nuclei. To avoid the effect of 
the source size and leave only the isospin effect, IMF mul- 
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Fig. 4. The average IMF multiplicity (top panel) and its vari- 
ance (bottom panel) per nucleon versus E* for various nuclei. 

tiplicity is divided by Aq. In Fig. 4, we present the scaled 
average IMF multiplicity and its variance versus E* . In 
this figure, one may clearly see an effect of the isospin and 
the source size on fragment multiplicities. 



4 Largest fragments in fragment partitions 

In multifragmentation, the mass number of the largest 
fragment A max may be used as an " order parameter" since 
it is directly related to the number of fragments produced 
during disintegration. We have analyzed its behavior with 
respect to the excitation energies. In order to remove the 
effect of source size and leave only an isospin effect, we 
have scaled it with the mass number of the source Aq. We 
have plotted A max /Ao as a function of excitation energy in 
Fig. 5 (top panel). The values of these quantities decrease 
first rapidly with E* , and then, beyond some value (E* ~ 
10 MeV/nucleon), decrease rather slowly. This suggests 
that the vaporization process becomes dominant around 
this point (see also an analysis in Ref. rjHj)- Another strik- 
ing result of this behavior is that all curves almost coincide 
in the transition region as displayed on the top panel in 
Fig. 5 (i.e. a universality of the behavior of A max ). We 
have also calculated the average A max variances in the 
multifragmentation stage for all nuclei at the same ex- 
citation energy range. The average variance of A max /Ao 
values are presented in the bottom panel of Fig. 5. We ob- 
serve from this figure that the maximum variance values 
for all sources are seen at an excitation energy of 4 — 5 
MeV/nucleon, which is exactly corresponding to the re- 
gion of the transition from compound- like channels to the 
full multifragmentation. In the following we clarify why 
this kind of behavior takes place. 




E*(MeV/nucleon) 

Fig. 5. < A max > and its variance a(A max ) divided by Ao 
versus E* for Au 197 , 5*n 124 , La 124 and Kr 7S . 

5 Temperature and bimodality 

Temperature of fragments is an important ingredient for 
any description of nuclear multifragmentation. Besides sta- 
tistical approaches, a temperature can even be introduced 
in dynamical (coalescence-like) processes of fragment pro- 
duction [231 - In SMM, the temperature of fragments Tf 
in separate partitions is defined from their energy balance 
according to the canonical prescription: 

"10+^ = l.5(n-l)T f + Y,(m t + E*(T f )) + E coul , (2) 

i 

where n is multiplicity of hot fragments and 'gas' par- 
ticles, mo and rrii are ground masses of initial nucleus 
and the fragments, E* are internal excitation energies of 
fragments, and E cou i is their Coulomb interaction energy. 
Masses and internal excitations are found according to 
their free energies (see eq. (1) and Ref. [13] )• The term 
with (n — 1) takes into account the center of mass con- 
straint, which is important for finite systems |24|. In this 
approach we allow for temperature to fluctuate from par- 
tition to partition, and we define the temperature of the 
system as the average one for the generated partitions. 
The methods involving fluctuations of temperature have 
been used in thermodynamics since long ago j^HI ■ Our def- 
inition of temperature has another advantage: It is fully 
consistent with the temperature commonly used in nuclear 
physics for isolated nuclei, that makes a natural connec- 
tion with the physics of compound nucleus. 

Figs. 6a and 6b show the results of the average tem- 
perature and its variance with respect to the excitation 
energy E* for all nuclei under consideration. While one 
may see a plateau in the excitation energy range 3 — 7 
MeV/nucleon in the upper panel of this figure, the vari- 
ance of temperature exhibits a maximum value at about 
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Fig. 7. The values of Tmin and Tmin as a function of the critical 
temperature T c for various nuclei. 



energy range. However, temperatures of these groups of 
fragments are essentially different, since the binding en- 
ergy effect is different, and a disintegration into a large 
number of fragments takes more energy. The transition to 
the energy consuming group A max < Aq /3 occurs because 
of the phase space domination, despite of a smaller tem- 
perature of this group. We will call it 'bimodality' from 
now on (see also discussion for other models in (2H|)- It is 
an intrinsic feature of the phase space population in the 
SMM, however, it was only recently realized that it can 
be manifested in different branches of the caloric curve 
and in momenta of fragment distributions (see Ref. \7\). 
The bimodality provides an explanation why the caloric 
curve may have a plateau- like (or even a 'back-bending') 
behavior and large fluctuations in the transition region. 



6 Influence of the critical temperature 



5 — 6 MeV for standard SMM parameters, in agreement 
with the maximum fluctuations of A max . The plateau- like 
caloric curve was reported long ago (e.g. [T^ITH?] '). It was 
confirmed by experimental results |26| . and by rather so- 
phisticated calculations of Fermionic Molecular Dynam- 
ics |2Zl and Antisymmetrized Molecular Dynamics |28| . 
The reason for large fluctuations of the temperature and 
Amax, can be clear from Figs. 6c and 6d. In these figures 
we show relative mass probability and temperature values 
for two groups of fragment partitions, which are defined as 
Amax > 2^4o/3 (associated with compound-like channels) 
and Amax < Aq/3 (full multifragmentation) . As can be 
seen in Fig. 6c, the probability of observing the first group 
events A ma x > 2^4.0 /3 decreases rapidly in the excitation 
energy range 2 — 6 MeV/nucleon, and this probability for 
the second group A ma x < A$ /3 increases within the same 



We have analyzed how the critical exponent r depends on 
T c , similar to Ref. ^^l, and on isospin of the sources. Dur- 
ing multifragmentation nuclei demonstrate critical behav- 
ior at considerably lower temperatures (T* w 5 — 6 MeV) 
than the critical temperature of nuclear matter T c (see 
also [ZEIE]), because of Coulomb and finite size effects. 
As was shown in many papers (see e.g. f^H] and references 
in) the critical behavior can manifest as scaling of some ob- 
servables around the critical point T* . In the present work 
we are interested in information about the critical temper- 
ature T c , which is high enough to be extracted through the 
scaling in finite nuclei. However, this information can be 
obtained by looking at other features of fragment distribu- 
tions, since the critical temperature influences the surface 
tension of fragments. We have extracted the minimum val- 
ues of t versus T c and the corresponding temperatures of 
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the systems, as shown in Fig. 7. The similar analysis of ex- 
perimental data was carried out by Karnaukhov et al. |17j 
for residues produced after intranuclear cascade in Au tar- 
get. One can see from Fig. 7 that both the r m i n and T m j„ 
the temperature corresponding to this minimum, increase 
with T c . However, nuclei with lower isospin such as La 
and Kr 7s exhibit much lower r m i„ values. It is because 
of the proton-rich nuclei are less stable, and they disap- 
pear rapidly with the excitation energy, this corresponds 
to a less pronounced U-shape mass distribution. This ef- 
fect may be important for explanation of some data, since 
sources with different isospin may manifest different T m j„ . 
The decreasing T m i n for lower T c is explained by a very 
fast decreasing of the surface tension with temperature 
in this case. Nevertheless, all parameters increase with T c , 
and tend to saturate at T c — ► oo, corresponding to the case 
of the temperature-independent surface. In this case only 
the translational and bulk entropies of fragments, but not 
the surface entropy, influence the probability of fragment 
formation. 



7 Influence of the symmetry energy 

Now let us turn to the analysis of the symmetry energy 
of fragments, and its influence on fragment partitions. As 
we discussed, the symmetry energy for fragments is de- 
fined in SMM as E s y ,n = 7 (A - 2Z) 2 /A, where 7 is a 
phenomenological coefficient. As an initial approximation 
for the SMM calculations we assume 7=25 MeV, corre- 
sponding to the mass formula of cold nuclei. We have also 
performed the calculations for 7=8 MeV for Sn 124 and 
La 124 sources. In Fig. 8 we demonstrate the caloric curve, 
the mean maximum mass of fragments in partitions, and r 
parameters for different assumptions on the symmetry en- 
ergy. The results for 7=8 MeV are only slightly different 
from those obtained for the standard SMM assumption 
7=25 MeV. From this figure, one can see a small decrease 
of temperature caused by the involvement of the fragments 
with unusual neutron and proton numbers into the dom- 
inating partitions. Therefore, the average characteristics 
of produced hot fragments are not very sensitive to the 
symmetry energy, and special efforts should be taken in 
order to single out this effect. 

The main effect of the symmetry energy is manifested 
in charge and mass variances of the hot fragments. In Fig. 

9 we have shown the relative mass distributions of primary 
hot (top panels) and final cold fragments (after secondary 
deexcitation, middle and bottom panels) with Z — 5 and 

10 at an excitation energy of 5 MeV/n for 7 = 25 and 8 
MeV. The distribution of primary hot fragments becomes 
rather broad for 7 = 8 MeV and its maximum lies on the 
neutron-rich side. 



8 Evaporation of hot fragments with isospin 
effects 

The primary fragments arc supposed to decay by sec- 
ondary evaporation of light particles, or by the Fermi- 
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Fig. 8. Variation of the caloric curves (top panel), A ma x/Ao 
(middle) and r (bottom) of the nuclei Sn 124 and La 124 with 
symmetry energy versus excitation energy. 

break-up |20U13| . However, these codes use standard mass 
formulae obtained from fitting masses of cold isolated nu- 
clei. If hot fragments in the freeze-out have smaller values 
of 7, their masses in the beginning of the secondary deex- 
citation will be different, and this effect should be taken 
into account in the evaporation process. Sequential evap- 
oration is considered only for large nuclei (A >16) evap- 
orating lightest particles (n,p,d,t, 3 Ke,a). We believe that 
we can estimate the effect of the symmetry energy evo- 
lution during the sequential evaporation by including the 
following prescription. There are two different regimes in 
this process. First, in the case when the internal excitation 
energy of this nucleus is large enough (E* /A > lMeV), we 
take for hot fragments the standard liquid drop masses mid 
adopted in the SMM, as follows 

m AZ = midi'j) = m n N + m p Z - AW + 7 , . 

,D 42/3 , 3e 2 Z 2 W 

where m„ and m p are masses of free neutron and proton. 
In the second regime corresponding to the lower excita- 
tion energies, we adopt a smooth transition to standard 
experimental masses with shell effects (m s t, taken from 
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Fig. 9. fsotopic distribution of primary and secondary frag- 
ments for Z — 5 and Z = 10 versus mass number A at an ex- 
citation energy of 5 MeV/nucleon for Au 197 , Sn 124 and La 124 . 
At 7=8 MeV the results after standard secondary decay are 
shown with dotted lines, and dashed lines denote the results 
of new evaporation depending on the symmetry energy(see the 
text). 



to the standard code [23j and the second one is with the 
above described version taking into account the symmetry 
energy (mass) evolution during evaporation. The standard 
deexcitation leads to narrowing the distributions and to 
concentration of isotopes closer to the 0— stability line by 
making the final distributions of the case of 7=8 MeV sim- 
ilar to the case of 7=25 MeV. However, a sensitivity to the 
initial values of 7 remains, as one can see from distribu- 
tions of cold fragments with Z=5 and 10. This difference 
in distributions is much more pronounced if we use the 
new evaporation. The final isotope distributions are con- 
siderably wider, and they are shifted to the neutron rich 
side, i.e. the produced fragments are neutron rich. This ef- 
fect has a simple explanation: By using the experimental 
masses at all steps of evaporation we suppress emission 
of charged particles by both the binding energy and the 
Coulomb barrier. Whereas, in the case of small 7 in the 
beginning of evaporation, the binding energy essentially 
favors an emission of charged particles. When the nucleus 
is cooled sufficiently down to restore the normal symmetry 
energy, the remaining excitation is rather low (E* /A <1 
MeV) for the nucleus to evaporate many neutrons. 

The difference between the two kind of evaporation 
calculations gives us a measure of uncertainty expected 
presently in the model. This uncertainty can be dimin- 
ished by comparison with experiments. It is encouraging, 
however, that all final isotope distributions have some de- 
pendence on the 7 parameter. This sensitivity gives us 
a chance to estimate the symmetry energy by looking at 
isotope characteristics. As was shown earlier the actual 
7 parameter of the symmetry energy can be determined 
via an isoscaling analysis [1231 • An analysis of N/Z ratio 
of produced fragments can also be used for this purpose. 
Recently one may see some experimental evidences for es- 
sential decreasing of the symmetry energy of fragments 
with temperature in multifragmentation |31II32| . The con- 
sequences of this decreasing are very important for astro- 
physical processes [2]. 



nuclear tables) with the following linear dependence, 

mAZ = mid(l) • x + m st ■ (1 - x), (4) 

where x = (3E*/A (J3= lMeV" 1 ) and x <1. The excitation 
energy E* is always determined from the energy balance 
taking into account the mass thaz at the given excitation. 
This mass correction was included in a new evaporation 
code developed on the basis of the old model [2U], tak- 
ing into account the conservation of energy, momentum, 
mass and charge number. We have checked that the new 
evaporation at 7=25 MeV leads to the results very close 
to those of the standard evaporation. 

Generally, the secondary deexcitation pushes the iso- 
topes towards the value of stability, however, the final 
distributions depend on the initial distributions for hot 
fragments. One can see that the results can also depend 
on whether the symmetry energy evolves during the evap- 
oration or not. In Fig. 9 we demonstrate the results of 
two kind of calculations. First one is carried out according 



9 Conclusions 

In summary, we have shown the effects of isospin, crit- 
ical temperature and symmetry energy on the fragment 
distribution in multifragmentation. For this purpose we 
have analyzed different nuclei with various neutron-to- 
proton ratio on the basis of SMM. The critical temper- 
ature of nuclear matter influences the fragment produc- 
tion in multifragmentation of nuclei through the surface 
energy, while the symmetry energy determines directly 
the neutron richness of the produced fragments. Effects 
of the critical temperature can be observed in the power 
law fitting parametrization of the fragment yields by find- 
ing r parameter. We have also shown the effect of the 
isospin and neutron excess of sources on the fragment 
distributions and on the r parametrization. By selecting 
partitions according to the maximum fragment charge we 
have demonstrated a bimodality as an essential feature of 
the phase transition in finite nuclear systems. This fea- 
ture may allow for identification of this phase transition 
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with the first order one. It has been found out that the 
symmetry energy of the hot fragments produced in the 
statistical freeze-out is very important for isotope distri- 
butions, but its influence is not very large on the mean 
fragment mass distributions after multifragmentation. We 
have shown that the symmetry energy effect on isotope 
distributions can survive after secondary deexcitation. An 
extraction of this symmetry energy from the data is im- 
portant for nuclear astrophysical studies. 
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